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ABSTRACT 

A short gamma-ray burst GRB 090510 detected by Fermi shows an extra 
spectral component between 10 MeV and 30 GeV, an addition to a more usual 
low-energy (< 10 MeV) Band component. In general, such an extra component 
could originate from accelerated protons. In particular, inverse Compton emis- 
sion from secondary electron-positron pairs and proton synchrotron emission are 
competitive models for reproducing the hard spectrum of the extra component in 
GRB 090510. Here, using Monte Carlo simulations, we test the hadronic scenar- 
ios against the observed properties. To reproduce the extra component around 
GeV with these models, the proton injection isotropic-equivalent luminosity is re- 
quired to be larger than 10^^ erg/ s. Such large proton luminosities are a challenge 
for the hadronic models. 

Subject headings: cosmic rays — gamma rays: bursts — gamma rays: theory — 
radiation mechanisms: nonthermal 

1. Introduction 

The very high bulk Lorentz factor and high magnetic field strengths in outflows of 
jamma-ray burst (GRB) make GRBs potential sources of ultra-high-energy cosmic rays 
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fUHECRs: IWaxman lllQQSl : IVietri lllQQSl : ICialis fc Pelletier Il2004f ). Recent time -resolved de- 



tection s of GeV photons with Fermi/LAT from GRB 080916C and GRB 090510 flAbdo et al 



200930) reveal that the minimum Lorentz factor Fmiru required to make the sources op- 
tically thin for such GeV photons, should be > 1000. While those high values of Fmin 
are favorable to accelerate protons to very high energies, it makes difficult to produce the 
photomeso n- induced secondary photons and neutrinos frequently discussed by many au- 



thors (e.g-.lRachen fc Meszaros 



Asano fc Nagataki 



2006 



1998 



Dernier 



2002; Asano 



20051 : iDermer fc Atoyan 



2006 



Gupta fc Zhang 1 120071 ). This is because, given the variability 



timescale At, photon luminosity, and spectrum, a higher bulk Lorentz factor F decreases 
the photon number density as oc F~^. The wide-b and spectra (10 ke V - GeV) of GRB 



080916C are well fitted by a smooth Band function (IBand et al. Ill993l ). and in that case 



the lack of evidence for an extra spectral component is consistent with a low efficiency of 
photomeson production due to high F. 

However, from 0.5 s to 1.0 s after the trigger time (Tq), the short GRB 090510 does 
exhibit a very significant (> 5a) spectral deviation to the standard Band function. With the 
parameters of the standard Band function at £peak = 3.94 ±0.28 MeV, a = —0.58 ±0.06, and 
13 = —2.83 ± 0.20, the excess is adequately fit with an additional power law of index —1.6. 
The power law extends from the lowest energy in Fermi/ GBM (i.e., 8 keV), to at least 31 
GeV in Fermi/LAT. The extra-component counts for about 37% of the total fiuence and the 
powerlaw dominates the standard integrated-Band spectrum up to a few tens of keV and 
above 10 MeV. This additional component is the most intense from To ± 0.6 s to To ± 0.8 
s. There is no clear evidence of existence of an extra component at times after Tq ± 0.8 s, 
however there is not enough statistical significance in the data to fully reject this hypothesis. 
The onset of the main emission above 100 MeV is delayed about 200 ms compared to the 
main spike from 8 keV up to few MeV. No lags have been found below 1 MeV with respect to 
the lowest energy band 8-40keV over the whole light curve ( To to Tq ± 1.5 s). Above 1 MeV, 
lags increase progressively to reach 248 ms and remaining constant after 40 MeV. This extra 
comp onent could be of le ptonic origin, e.g., external Compton emission from internal shocks 
(e.g., iToma et al. I l2009l ) or synchrotron self -Compton emission fro m the reverse shock or 



forward shock in the early afterglow phase (Ghirlanda et al. II2009I ) as discussed for GRB 
941017 JGranot fc Guetta Il2003l : 1pe'er fc Waxman Il2004h . 



In this Letter we consider an alternative possibility, nam ely hadronic models, in par- 



ticular photomeson c ascade and proton synchrotron models (jVietri 1 119971 : iTotani I Il998 



Razzaque et al. II2009I ). for representing the extra spectral component of GRB 090510, which 
is present only in the prompt emission between To and To ± 0.8 s. Although the cascade pro- 
cesses initiated by p ± 7 p/n + vr^/vr"'" are complicated, the resultant photon signatures of 
proton cascades mostly appear as synchrotron or inverse Compton emission from secondary 
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electron-positron pairs produced via 7 + 7 — e"*" + e . In such pair cascade proce sses the 



effective injection index of secondary pairs tends to be about —2 (e.g. ICoppi Ill992l ) so that 



the synchrotron radiation from secondary pairs yields a fiat spectrum in an e'^N{e) plot, 



inverse Compton component from secondary pairs can harden the spect rum as lAsanq et al. 



rum as 


Asano et al. 


Asano et al. 


( 


2009) 



(see also lAsano fc Inoue 1120071 ). we constrain the hadronic models in this Letter. 

After a short discussion of our methods in |J2], we discuss the possibility of photon 
emission due to accelerated protons for GRB 090510 in ^ and summarize the results in ^ 



2. Model and Methods 



In order to constrain hadronic models for the observed extra spectral component of 
GRB 090510, we simulate the photon emission with photomeson-induced cascade processes 
from ul trarelativistic outflo ws of bulk Lorentz factor F. The numerical code i s the same as 
t hat inlAsano et al. I (120091). which was niatured via a series of GRB studies ( lAsano 112005 



Asano &: Nagataki 



2006 



Asano fc Inoue 1120071 ). so that we omit the detailed explanation 



for the code. As an optimistic case, which would lead to hadronic cascades, we assume that 
protons are injected with a power law energy distribution oc i^'^exp {—E^/E^^.^^ above 10 
GeV in the outflow rest frame at radii R from the central engine. The acceleration timescale 
of a proton is parameterized as tacc = ^Rl/c, where i?L is the Larmor radius of the proton. 
The maximum proton energy -Emax is determined by equating the acceleration timescale and 
the expansion timescale of the outflow texp = -R/(cr) or the cooling timescale due to proton 
synchrotron or photomeson production. 

We assume that the main low-energy component, which is fitted by a Band function, 
is of leptonic origin lAbdo et al. I (l2009bl ). which we do not discuss further. The photon 
energy distribution in the outflow rest frame is estimated from the time-averaged Band 
parameters from Tq + 0.5 s to Tq -|- 1.0 s (^peak = 3.9 MeV, a = —0.58, /3 = —2.83, and 
the isotropic-equivalent luminosity = 2.5 x 10^^ erg/s aX z = 0.903) and the numerically 
obtained extra component originating from protons. The energy distributions of photons 
and particles are simulated iteratively until they converge to a self-consistent spectrum. On 
a timescale comparable to texp, the photon density and the magnetic fields, etc., can be 
approximately taken as constant, and we can neglect the emission from particles beyond this 
timescale after their injection, because the photon density and the magnetic field decline. The 
physical processes taken into account are 1) photon emission processes of synchrotron and 
Klein-Nishina regime Compton scattering for electrons/positrons, protons, pions, and muons. 
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2) synchrotron self-absorption for electrons/positrons, 3) 77 pair production, 4) photomeson 
production from protons and neutrons, 5) Bethe-Heitler pair production (p+7 p+e~^+e~), 
and 6) decays of pions, and muons. The code can also output the spectra of neutrinos from 
pions and muons. 

There are five model parameters: the proton acceleration parameter ^, the proton in- 
jection radius R, the bulk Lorentz factor F, the energy density of the magnetic field Ub, and 
the injection luminosity of accelerated protons Lp. The last two parameters are normalized 
by the photon energy density of the Band component = L^/(47rc-R^r^) or itself as 
Ub/U^ and Lp/L^. 



3. Hadronic Emission 

Within the time interval between Tq + 0.5 s and Tq + 1.0 s, two high-energy photons 
of 3.4 GeV and 31 GeV were detected. These photons belong to the extra component of 
the prompt emission, whose spectrum is well fitted by a power law of —1.6, on which we 
concentrate here. (There is in addition a tail extending to ~ 200 s, attributable to an 
afterglow, which we do not discuss here.) Since the photon statistics above GeV are not 
enough, we do not attempt to carry out a detailed fitting of the spectrum. Alternatively, 
we search for parameter sets that reproduce a comparable fiuence of the extra component 
to that of the Band component. 

First, we consider photomeson cascades, for the most extreme case ^ = 1. Fig. [1] shows 
two examples of photon spectra, in which the flux of the hadronic component at 3.4 GeV 
is comparable to the Band component. The model spectrum with Ub/U^ = 10~^ seems 
consistent with the observed spectrum. The second peak at ~ GeV is mainly due to 77 pair 
absorption. The assumed value of Lp/L^ is 200, which is quite large, so that the proton 
injection luminosity should be larger than 5 x 10^^ erg/s in this case. On the other hand, if 
we adopt a stronger magnetic field such as Ub/U^ = 10~\ the required amount of protons 
can be suppressed to a lower value of Lp/L^ = 30 (Lp ~ 7 x 10^^ erg/s). However, the 
dominance of synchrotron radiation from secondary pairs makes in this case the spectrum 
too soft below about 100 keV. Such a large deviation from the Band function around 10 keV 
is not seen in GRB 090510. 

The various timescales for the case of Ub/U^ = 10^^ are plotted in Fig. [2J Apparently, 
the maximum energy of protons -Emax is determined by the condition tgxp = ^acc- The low 
efficiency of hadronic cascade is attributed to the much longer cooling timescale due to 
photopion production t.„ than the expansion timescale of the outflow texp even at Ep = E^^slx- 
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The target photons for pion production below ~ 3 x 10^^ eV are the Band component. 
Above this energy, protons interact with UV-photons below the observed energy range of 
Fermi, where secondary synchrotron photons dominate. In our simulation the synchrotron 
self-absorption makes a spectral peak of eN{e) plot at ~ 40 eV, which corresponds to the 
typical target photon energy for protons of ~ 10^^ eV, where becomes minimum. In spite 
of the enhancement of the pion production efficiency due to secondary photons, we still need 
a large proton luminosity. The curves in Fig. [2] can be shifted by changes of R, F, and 
Ub/U^- The increase of -Emax due to the decrease of tacc oc B^^ with increasing magnetic 
fields leads to a higher efficiency of photomeson production. As shown in Fig. [1], however, 
we need to decrease the magnetic field to values as low as Ub/U^ = 10"^ in order to harden 
the spectrum, since the Klein-Nishina effect is crucial. 

Let us consider protons that interact with photons of ^peak- The cooling timescale for 
such protons behaves as t-,, oc n~^, where the photon density = U^/ (^pcak/r)- If we adopt 
a smaller value of F, the pion production efficiency would be improved as texp/^7r -R^^F^^. 
However, we should note that there is a lower limit to F, which is required to make the 
source optically thin to GeV photons. Given the photon luminosity and spectral shape, this 
minimum Lorentz fa ctor can be estimated as shown in the online supporting materials in 
Abdo et al. ] J2009al lb[). The adopted value F = 1500 in Fig. [1] is close to the lower limit for 



R = 10^^ cm. Even if we take a larger R, the lower limit of the Lorentz factor oc R^/^l^~^^ 
does not decrease drastically (since /3 ~ —3, Fmin oc Although a smaller R <C 10^^ 

cm and a slightly larger F can reduce the required amount of protons, the typical variability 
timescale i?/(cF^) in Fig. [H (i? = 10^^ cm and F = 1500) is already as small as 1.5 ms. The 
poor photon statistics in GRB 090510 cannot constrain well the variability timescale, but 
the required proton energy may remain very large even for i?/(cF^) <^ 1 ms (note that in 



Abdo et al. I (j2009bl ) the variability timescales are estimated as ~ 10 ms, which are adopted 



to obtain the minimum Lorentz factor). 

An alternative hadronic scenario is a proton synchrotron model with very strong mag- 
netic field. In this model, we need to avoid the pair cascades initiated by photomeson produc- 
tion, because the synchrotron spectrum from secondary pairs is too soft as discussed above. 
Therefore, we adopt ^ = 10^ to suppress -Emax- Examples for this model in Fig. [3] require the 
same order of magnitude ratios of Lp/L^ as the pair cascade model. Considering the typical 
energy of protons emitting GeV photons, oc (i?F)~^/^, the efficiency of proton synchrotron 
can be enhanced by a decrease of F, according to texpAsyn oc T^^R^^/'^(Ub/UjY^^. However, 
even for F = 3000 in Fig. [31 the pair cascades due to photomeson production cannot be 
neglected in the low-energy regions. The photon spectrum below 100 keV becomes too soft 
again for F = 2000 owing to the pair cascade. If we can neglect the observed variability 
timescale < 10 ms, it seems that a much larger R and slightly smaller F could decrease the 
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Fig. 1. — Photon spectra (bold curve) for ^ = 1, R = 10^^ cm, and F = 1500. The red 
curve is for Ub/U^ = 10~^ and Lp/L^ = 200, and the black curve is for Ub/U^ = 10~^ 
and Lp/L^ = 30. Fine red curves denote separately pair synchrotron and inverse Compton 
without the absorption effects. The fine dashed curve is the Band component. The vertical 
blue line denotes the energy of the 3.4 GeV photon. 

required luminosity of protons. Under the conservative assumption of f/^/f/p < 1, however, 
a decrease of Up leads to a low efficiency of proton synchrotron due to the low Ub- Consid- 
ering the dependence of texp/4yn oc T^^R^^^'^{UB/Uy)'^^^, we cannot avoid the requirement 
of a large proton luminosity even in the proton synchrotron scenario. 

In the above, we have sought to reproduce the flux at 3.4 GeV by photons originating 
from protons. To extend the extra component as far as 31 GeV, a much larger F is required, 
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t [s] R=10^4cm, r=1500, Ub/U^=10-^ Lp/L^=200 




Ep [eV] 

Fig. 2. — Various timescales for the model of Ub/U^ = 10~^ in Fig. [1) proton cooling 
timescales due to photomeson production {t.,,), Bethe-Heitler pair production (tBH), syn- 
chrotron emission (tsyn), and inverse Compton emission (tic), proton acceleration timescale 
(tacc)) and expansion timescale of the outflow (texp) in the outflow rest frame. 

which would decrease the efficiency of photopion production, (for example, L^/L^ = 1000 is 
required to obtain sufficient flux at 31 GeV for ^ = 1, F = 3000, Ub/Uj = 10"^, and R = 10^^ 
cm). Thus, we need an isotropic equivalent proton luminosity in excess of 10^^ erg/s in order 
to explain the photon flux and hard spectrum of the extra component, independently of the 
model details. 



Another aspect of this burst is the reported ~ 0.1-0.2 s delayed onset of the LAT high- 
energy photons relative t o the GBM trigger. This delayed onset is a common property of 
other Fermi/LAT GRBs ( lAbdo et al. Il2009al etc.). Within the hadronic models, this delay 
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Fig. 3. — Photon spectra (bold curve) for ^ = 10^ and R = 10^^ cm. The red curve 
is for r = 3000, Ub/U^ = 300, and Lp/L^ = 300, and the black curve is for T = 2000, 
Ub/U^ = 200, and Lp/L^ = 100. Fine red curves denote separately proton synchrotron, pair 
synchrotron, and muon synchrotron without the absorption effects for F = 2000. The fine 
dashed curve is the Band component. 



could in principle be attributed to the acceleration timescale of protons, which is limited by 
the shell expansion timescale. This would suggest that the timescale of the delayed onset 
is comparable to the variability timescale, which may be comparable to or shorter than 10 
ms. However, since the cooling timescale is much longer than the shell expansion timescale 
in our models (which is equivalent to the requirement Lp/L^ ^ 1), we cannot deduce 
the delay timescale in this manner from the models considered here. In order to explain 
the delayed onset, we would need additional assumptions, an example of which might be 
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spatially different origins of the Band and the extra component. However, with only two 
photons above 1 GeV, it is premature to engage here in such additional modeling. 

In Fig. m we have also plotted the neutrino spectra (sum over all species) obtained 
from the model for Ub/U^ = 10~^ in Fig. [1] (photopion- induced pair cascade model) and 
for r = 3000 in Fig. [3] (proton synchrotron model). For the pair cascade case, the very 
large proton luminosity leads to a neutrino luminosity comparable to the photon luminosity. 
The neutrino number flux is estimated as ~ 6 km~^ s~^ for each energy decade from 10^^ 



to 10^^ eV. Considering the effective area of the IceCube detector (lAbbasi et al. II2009I ). we 
may not expect neutrino detections from this burst (the upper limit of the muon neutrino 
fiuence for the "naked-eye" GRB 080319B by the IceCube detector is ~ 10~^erg cm~^ for 
< 2.2 X 10^^ eV). For the proton synchrotron model, the strong magnetic field efficiently 
cools pions and muons via synchrotron radiation so that the highest neutrino energies are 
suppressed, compared to the pair cascade model. 



4. Conclusions 



The detection of GeV photons from the short GRB 090510 requires a very large Lorentz 
factor, which leads to a low efficiency of photomeson production. Adopting the observed 
spectrum, we have simulated the photon emission due to accelerated protons to produce the 
extra spectral component around GeV. We have shown that this mechanism is capable of 
reproducing the flux in this extra component, during the approximate time interval between 
the trigger and 0.8 s when this component is present. However, because of the low efficiency 
of photopion production, a very large proton isotropic-equivalent luminosity of > 10^^ erg/s 
is required to produce GeV photons from electron-positron pairs. While we have assumed 
a power-law proton spectrum with an exponential cutoff for simplicity, the energy range, in 
which photopion production predominantly contributes to the secondary photon emissions, 
is narrow as shown in Fig. [2l The required proton luminosity we suggest here practically 
represents the required normalization of the proton flux at this energy range (~ 10^^ eV 
for the example in Fig. [2]). Even if we change the proton injection spectrum, the effect 
on the photon spectrum would be negligible as long as the proton amount in that energy 
range is the same. Fixing the proton amount in this energy range, of course, steeper proton 
spectra enhance the total injection luminosity. In principle, a lower proton luminosity as- 
sociated with a stronger magnetic field could reproduce the flux around GeV, provided the 
spectral modifications below 100 keV due to secondary pairs are weak. This remains to be 
investigated, but in the absence of some physical mechanism to suppress the emission of low- 
energy photons from secondary pairs, the spectrum predicted from such a strong field model 
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Fig. 4. — Neutrino spectra (black bold curves) for the pair cascade model (F = 1500, 
,^ = 1, Lp/L^ = 200, Ub/U^ = 10~^, and R = 10^'^ cm) and the proton synchrotron model 
(F = 3000, ^ = 10^ Lp/L^ = 300, Ub/U^ = 300, and R = 10^^ cm). Fine curves denote 
separately neutrinos from pion decay and muon decay. 



would contradict the observations. Also for the proton synchrotron model with UB/U-y ^ 1 
discussed above, the required proton isotropic-equivalent luminosity is > 10^^ erg/s. 



Such "proton-dominated" GRB models (e.g.. lAsano et al. 1120091 ) are favorable for sce- 
narios of ultra-high-energy cosmic-ray generation by GRBs, but as in this case too, they 
require very high isotropic-equivalent proton energies > 10^^ erg/s, which are challenging for 
such hadronic models. However, we should note that even the photon isotropic equivalent 
luminosities of quite a few bursts have been found to be of that order of magnitude (e.g., 
GRB 080913, etc.), so that the hadronic interpretation of the extra GeV component does 
not exceed such observed photon luminosities in other bursts. Jet coUimation effects can sig- 
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nificantly alleviate such energetic problems. Unfortunately, in the case of GRB 090510 the 
afterglow is brief, and there is so far no evidence for or against any jet collimation. Adopting 
1/500 > 1/r as the collimation factor, the required energy may strain a conventional NS-BH 
merger model, but is not huge. The strain would be less on a core collapse model, but this 
model has not solved the short- duration issue yet. Therefore, the extreme energetics of such 
hadronic outflow models remain a substantial challenge. 

In order to distinguish models of in ternal shock origin including hadronic models and 
afterglow origin (iGhirlanda et al. I l2009l ) for GeV photons, better photon statistics above 
GeV are desirable. Although the correlation function for this burst shows a clear correl ation 
between GEM (250 keV-3MeV) and LAT (> 0.1 GeV) light curves JAbdo et al. Il2009bh . the 
typical photon number above GeV in a pulse is a few. Closer GRBs than GRB 090510 at 
least by a factor of 3, from which we expect more than 10 photons in a pulse, may strengthen 
the correlation, but we should note that the GRB rate may decrease with increasing photon 
statistics (the present detection rate with Fermi LAT is ~ 10 GRBs yr^^). 



We thank the Fermi GBM/LAT collaborations, especially S. Razzaque, C. Dermer, and 
J. Granot, for valuable discussions, and acknowledge partial support from NSF PHY-0757155 
and NASA NNX 08AL40G (P.M.). Finally, we appreciate the anonymous referee. 
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